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INTRODUCTION 


Many  naval  systems  use  time  difference  methods  to  estimate 
position.  While  the  arrival  time  of  a transmitted  signal  de- 
pends on  the  transmission  time  and  the  length  of  the  trans- 
mission path,  differences  of  arrival  times  do  not  depend  on  the 
unknown  time  at  which  the  transmission  takes  place.  Conse- 
quently, estimates  of  the  position  of  a transmitter  or  receiver 
are  based  on  observed  time  differences,  and  the  quality  of  the 
estimate  depends  on  the  errors  in  these  quantities,  that  is,  on 
the  joint  time  difference  error  density. 

Although  a knowledge  of  this  density  function  is  required 
for  a quantitative  treatment  of  any  aspect  of  the  estimation 
problem,  it  will  seldom  be  directly  available.  When  the  system 
is  calibrated,  measurements  performed  on  transmissions  between 
known  locations  or  over  known  paths  provide  data  on  the  errors 
in  the  observed  arrival  times.  This  process  yields  a joint 
density  for  the  errors  in  the  arrival  times,  and  the  corres- 
ponding density  for  the  differences  in  the  arrival  times  must  be 
derived  from  it.  A formulation  of  the  density  transformation 
and  some  examples  of  its  application  are  given  in  this  research 
contribution . 


FORMULATION  OF  THE  TRANSFORMATION 


If  there  are  n observed  arrival  times  there  are  n(n-l) 
time  differences.  Of  these,  only  (n-1)  may  be  chosen  indepen- 
dently. It  will  therefore  be  assumed  that  one  time  is  selected 
as  a reference  and  that  the  remaining  times,  measured  with  re- 
spect to  this  reference,  constitute  the  (n-1)  differences  of 
interest . 


Since  taking  differences  is  a linear  operation  which  in- 
troduces correlation,  matrix  notation  is  convenient,  particu- 
larly if  related  computer  programs  are  written  in  a language 
suitable  for  processing  numerical  arrays.  We  therefore  let 


denote  the  observed  and  true  values  of  the  arrival  times,  re- 
spectively. The  observed  and  true  values  of  the  differences  in 
arrival  times,  where  the  kttn  time  is  selected  as  reference,  are  then 
given  by  and  , respectively,  where  A^  is  the  (n-1)  x n 

matrix  derived  from  the  n x n identity  matrix  by  replacing  each 
zero  element  of  the  kfc^  column  by  -1  and  deleting  the  kfc^  row. 

Thus,  for  example,  if  there  are  five  observations  and  the  fourth 
time  is  selected  as  reference,  the  observed  time  differences  are 
given  by 


— 

— 

T1 

" T4 

i 

0 

0 

-1 

0 

T2 

- T4 

0 

1 

0 

-1 

0 

T3 

" T4 

0 

0 

1 

-1 

0 

T5 

1 

rr 

E- 

1 

0 

0 

0 

-1 

‘J 

Ll 


1 4 


A41 


T 


5 


-2- 


Representing  arrival  time  errors  by 


[*i 


L 2 ' 


n r 


= T - t , 

and  the  time  difference  errors  with  the  kth  time  as  reference  by 


-k  ‘ [£ 


lk  e2k  e (k-1 ) k e (k+1 ) k ’**  enk 


]' 


for  k = 1,  2,...,  n,  we  have 

£k  = Ak!  - Akt 


-k-  * 


(1) 


This  is  the  required  error  transformation.  It  permits  the 


those 

of  the  arrival 

E(ek) 

* -k- 

where 

U = 

["l  1 

“i  = 

E(ti) 

and 

Ee  = 

- ek 

-k 

where 

E . and 

u — ck 

error 

vectors  £ 

and 

are  given 

by 

(2) 


iBi > 2 , • • • f k f • . • f n 


(3) 


are  the  covariance  matrices  for  the 


cov ( e . , c . ) = o . 
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the  above  expressions  may  be  rewritten: 


E(eikl  = *.± 


(4) 


and 


^^ik^jk*  = °ij  “ °ik 


kj  + " kk 


(5) 


i # j ~ , • • « , ( k-1 ) , ( k + 1 ) / 1 • • i n * 

If  arrival  time  errors  are  uncorrelated,  then 


cov (e . . ,e . , ) 
lk  jk 


9 9 

*-  .dm 

a ■ 6 . . +o. 

l l j k 


(6) 


The  subsequent  treatment  of  the  case  where  the  arrival  time 
errors  are  jointly  normal  will  make  use  of  these  results. 


In  formulating  the  joint  error  density  transformat  ion  for 
the  general  case,  we  note  that  although  the  mapping  £ * is 


not  one-to-one,  the  mapping  t-  -*  e+  is,  where  G+  is  defined  by 

k — k 


Sk  ' [e 


p P 

lk  2k 


(k-l)k  : k e (k+1) k **'  °nk 


]' 


The  inverse  mapping 


^k 


t has  the  form 


el  °lk  + ;k 


e2  e2k  + ^k 


ek  = 'k 


c = e . + e . 
n nk  k 
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r 


and  the  Jacobian  ot  the  transformation  is  unity.  Therefore,  the 
required  transformation  is 


Vs*’ 


e.  ulk  e2k  (k-l)k  ik+1 ) k °nk 
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Pe*lelk  e2k 


• C i . . . VI  6.1  . i \ >,  •••??•)  il  u 

(k-1) k ( k ^ 1 ^ k nk 
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where 
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l 2 * " ! k 
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Ik- 1 ) k 
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t ) 

n 


is  the  joint  error  density  for  the  arrival  times 
If  these  errors  are  independent,  then 


pe  <ek> 

-k  K 


liO  n , P le . . +u)  du 

, i 1 : . l k 

k i*k 


IS) 


These  expressions  provide  the  formal  transformations  for 
obtaining  the  joint  time  difference  error  density  where  the 
joint  error  density  for  the  arrival  times  has  a known  mathe- 
matical form. 


REFERENCE  TIME  SELECTION 


Min-  thl!  section  we  consider  the  effect  of  changing  the 
nme  selected  ,s  reference.  For  two  such  selections? 


£rit£ 


e = A e 
-r  — r 1 


has  the  property 


-rk  -k  ~ Ar 


then 


That  is,  the  linear  operation  denoted  by 
differences  with  respect  to  the  rth  time 
with  respect  to  the  kth  time.  Since 


— rk  Provides  the  time 
from  those  measured 


A.  = T,  A = T,  T A 

-k  -kr  -r  -kr  -rk  -k  ' 


we  have 


-kr  -rk  -n-l 


where  I 


-n-l  ls  the  ) *(n-l ) identity  matrix,  or, 

Tkr  * • 

While  on.y  this  result  will  be  used  to  investigate  the  effect  of 
changing  the  reference  time,  we  note  that 


-k  — k ~ = 1 . + J 

K K -r  -r  -n-l  -n-l 
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is  an 


where  in-i  is  the  (n-l)x(n-l)  identify  matrix  and  Jn-1  is  an 

(n-l)x(n-l)  matrix  with  all  elements  equal  to  unity.  Since  it 
is  readily  shown  that 

(I  + J i )-1  = I . - i Jn  i , 

— n-l  — n-1  — n-l  n — n-1 

(A,A  ' )_1  = In  1 ~ “ J , 

— k— k -n-1  n —n-1 

exists  for  n > 0,  Trk  may  therefore  be  expressed  in  terms  of 
difference  operators  alone  by 

irk  * — r-k^  - — r— k 'Mk1'' 


where' 


-1  ' ' -1 

Ak  = Ak  (AkAk) 


is  the  generalized  inverse  of  the  (n-l)x  n matrix  Ak  . 

To  illustrate  the  effect  of  changing  the  reference  time, 
let  £*  represent  a sample  value  of  the  arrival  time  error.  If 
the  ktntime  is  selected  as  reference,  the  sample  values  for 
time  difference  errors  are  given  by 

ej  - 4 £*  . 

and  the  joint  error  probability  is  given  by  P (e*) . If, 

e,  k 

• t h ^ 

al ternati vely,  the  r time  is  selected  as  reference,  a 

different  time  difference  sample  results,  namely 

£*  = Ar  e*, 

and  the  joint  error  probability  is  then  given  by  PQ  (e*) . 

But,  since  e_r  = Trk  ek  and  since  the  absolute  value  of  the 
Jacobian  of  the  transformation  Trk  is  unity. 


NORMAL  ERRORS 


If  the  time  of  arrival  errors  are  jointly  normal,  that  is 

P (e)  ^ N(u,  L ) * 

then,  from  equations  2 and  3,  the  time  difference  error  density 
is  given  by 


’ek(-k)  v N(-k  -k 


A.  ) 
c — k 


and  the  elements  of  the  mean  vector  and  covariance  matrix  are 
given  by  equations  4 and  5,  respectively. 

If  the  arrival  time  errors  are  uncorrelated,  the  numerical 
inversion  of  the  covariance  matrix 

JL  ■ -k  4 

k 

may  be  avoided,  and  the  computation  of  error  probabilities  thus 
simplified.  For  this  case,  the  elements  of  £ are  given  by 
equation  6,  namely 

cov(eik'  ejk)  = °i  lSij  + °k  * 


S ince 


the  elements  of 


are  given  by 


(9) 


Therefore,  assuming  biases  have  been  removed  from  the  arrival 
time  errors,  the  required  time  difference  error  probabilities 
may  be  computed  by  means  of  the  expression 


K 

n- 


1 


exp  (-  j eR 


(10) 


where1 


K 


n-1 


<2tt) 


(n-l)/2 


-1/2 


= ( 2tt  ) 


(n-l)/2  1/2  T1 0 

i=l 


1Graybill,  F.A.,  "Introduction  to  Matrices  With  Applications  in 
Statistics,"  Wadsworth  Publishing  Co.,  Inc.,  Belmont,  Ca.,  1969, 
p.  184. 
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CONDITIONALLY  NORMAL  ERRORS 


In  general,  where  arrival  time  errors  are  non-normal, 
resort  must  be  made  to  the  transformations  of  equation  7 or  8 to 
obtain  the  required  time  difference  error  density,  A special 
case  of  practical  interest,  however,  may  be  formulated  in  terms 
of  densities  which  are  conditionally  normal.  It  may  be  possible 
to  assert  that  arrival  time  errors  are  normal  and  to  specify  the 
parameters  of  the  density  function,  subject  to  certain  condi- 
tions for  which  the  probabilities  of  occurrence  are  known.  An 
example  would  be  the  use  of  clocks  of  different  accuracies.  If 
the  error  introduced  by  each  clock  is  normal,  Am,  m *»  1,2,..., 
M',  denotes  the  event  that  clock  m is  used,  and  the  occurrence 
probabilities  P(Am)  are  known  for  each  m,  the  density  function 
for  a single  arrival  time  error  t ■ would  be 


21 

m=l 


M 


P(t 


A ) 
m 


P ( A ) 

m 


Clearly,  this  density  is  normal  only  if  the  use  of  one 
particular  clock  is  certain. 

For  the  general  case  we  may  postulate  Cm,  m-1,2,...,  M 

mutually  exclusive  and  exhaustive  conditions  for  which  P (i  |C  ) 

e — 1 m 

is  jointly  normal  and  may  be  denoted  by  N(0,  2J  i ) . The 
probabilities  of  occurrence  of  each  condition 


Pm  * P(C\J 
m ' m 


m= 1 , 2 , . . . , M 


will  be  assumed  known,  although,  in  practice,  estimates  of  these 
quantities  would  be  used. 

From  the  above  discussion  of  the  normal  case 

re„  ‘Sklcm>  - N(0.  1L  I 
-k  k t m 


where 
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V =“  A.  . A. 

Z-  -k  ■**■(  | m -k 

^k  | m 

and  the  required  joint  time  difference  error  density  is  given  by 


P„  <ei>  'v  U M P P (e.  1C  ) 
^k  ~k  m-1  m ?k  “k  m 


£ p„,  N|2'  , ) 
m km 


As  in  the  normal  case,  these  expressions  simplify  where 
the  arrival  times  are  uncorrelated  for  each  of  the  applicable 

conditions.  Elements  of  the  covariance  matrices  V V*  , 

— € m ' =e,  | 
km 


and 


-1 

k J m 

are 

then 

given  by 
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respectively.  Thus, 


» y*  A 

P_  (e.  ) = W'  P K | exp  (-  i — k — 

e.  -k  ^ m n-1  m v 2 e,  i 

-k  m km 


<1D 


where 


K ,1  =»  ( 2 TT  ) 

n-1  m ' 


- (n-1 ) /2  1/2 
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Care  must  be  taken  in  listing  the  mutually  exclusive 
conditions.  If,  in  the  above  example,  each  of  the  M'  clocks 
could  be  used  to  make  any  of  the  n time  measurements,  the 
total  number  of  possible  conditions  would  be  (M ’ ) l/(M ' -n ) ! . II 

there  are  M'  types  of  clock  and  more  than  n clocks  of  each 
type,  (M')n  conditions  would  be  possible. 

The  computation  of  time  difference  error  probabilities  may 
be  illustrated  by  some  examples.  Consider  the  case  of  three  in- 
dependent time  measurements  using  clocks  of  two  different  types, 
A and  B.  Each  type  of  clock  exhibits  normal  errors  with  zero 

0 7 

_ - j ■ . * j - *-  : 1 ..  mu « 


means  and  variances 


and  o! 


, respectively.  The  eight 


possible  cases  are  listed  below  together  with  occurrence 
probabilities  for  a simple  bimonial  model. 
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From  equation  11, 


(2k> 


L8 


exp 


1 
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which  for  this  case  may  be  written 


P (e  ) 
ek  -k 


e1  e1  e1  p rr+s  a-pA)3_(q+r+s)K2| 

q=0  r=0  s=0 


qrs 


exp  ( 


'-1 


i I 

'k  | qrs 


Implementation  of  this  formula  requires  only  the  computation  of 

y-i 

and  K . for  q,r,s  * 1,0.  Let  the  second 

ek|4rs  2|qrS 

observed  time  be  chosen  as  reference.  The  cases  denoted  above 
by  m»l  and  m*2  are  given  by  q*r»s=l  and  q«r*l,  s*0, 
respectively.  Thus, 
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required 
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Pl(1-p2)P, 

P1  U-P2>  (1-P3) 

(1~Pl,P2r3 
(1'r1)P2(l*PJ 
: 1-Pi ) (1'P2)p3 

1“pl) (1~p^) (1-p,) 
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V8 
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v 
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N(0, 


) 


V1 


may  be  written 

Pe  (£J  ' V1 


1 


where  the  covariance  matrices 


are  given  by 


qi'S 


V 

— * qts 


°iq  + °b(1~q) 


2r  * °b(1'r) 


°3S  * °b(1_s) 


for  q , r,  s ■ 1,  0.  The  diagonal  elements  are  therefore  e*'  , 

V ■*  - 1 

i * 1,2,3,  or  , and  the  values  of  ^ and  K „ 

b ^V,r.  ‘ ''rs 

may  be  computed  as  in  the  previous  example.  This  permits  the 
computation  of  ) for  each  case,  and,  hence,  the 

* k qrs 

value  of  r (e.  ) by  means  of  the  above  expression. 
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APPENDIX 


PROGRAM  NOTES  AND  LISTING 

The  number  of  time  observations  and  the  time  chosen  as 
reference  are  represented  by  means  of  the  binary  vector  KVEC. 
Thus,  for  three  observations  with  differences  taken  with  respect 
to  the  second 

KVEC  <-  1 0 1 . 


The  difference  operator  A,  is  then  given  by  the  function 
DELOP.  For  example, 

NVEC< 1 0 1 

DEL  OF' 

1 “1  0 

0 "1  1 


The  function  DEL  with  argument  H gives  the 

covariance  matrix  b 


Thus,  for  example. 


SIGEF' 

0.01  0 0 
0 0.01  0 
0 0 0.01 

DEL  SIGEF' 
0.02  0,01 
0.01  0.02 


A-l 


The  inverse  of  this  matrix  required  for  the  numerical 
evaluation  of  the  joint  density  function  is  provided  by  the 
function  INDEL.  Thus, 

INDEL  S 1 CEP 

66 . 66666667  "33 . 33333333 

3 3*3333  3 3 33  66*6 6666667 


This  function  is  therefore  equivalent  to  the  sequential 
application  of  the  defined  function  DEL  and  the  APL  primitive 


function 


in  many  systems, 
efficient . 


This  latter  function,  however,  is  unavailable 
and  in  any  case,  the  function  INDEL  is  more 


The  computation  of  time  difference  error  probabilities  for 
the  normal  case  is  by  means  of  the  function  MVNPROB  which  incor- 
porates all  of  the  above  functions.  The  left  argument  consists 
of  the  time  difference  errors  expressed  as  a vector.  The  right 
argument  may  consist  of  either  the  covariance  matrix  for  the 
time-of-arrival  errors  or  a vector  composed  of  the  diagonal 
(nonzero)  elements  of  the  matrix. 

The  function  CONPROB  effects  the  multiple  application  of 
MVNPROB  necessary  when  different  conditions  may  prevail.  The 
right  argument  of  this  function  consists  of  a single  matrix, 
each  row  of  which  is  equal  to  the  trace  of  the  covariance  matrix 
appropriate  to  one  condition.  The  weighting  of  the  outputs  of 
CONPROB  by  the  probability  of  occurrence  of  each  condition  is 
implemented  by  the  function  WTPROB.  It  follows  that  in  a 
situation  where  conditional  normality  is  applicable,  the  joint 
probability  of  an  arbitrary  vector  of  time  differences  may  be 
computed  by  means  of  the  single  entry 


PM  WTPROB  (X  CONTROB  TM) 

where  PM  is  a vector  of  occurrence  probabilities,  TM  is  the 
matrix  fabricated  row  by  row  from  the  traces  of  the  covariance 
matrices  applicable  to  each  condition,  and  X is  the  time  dif- 
ference error  vector  for  which  a probability  is  required. 
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VDEL0PCD3? 
v Z«-DELOP  i N 
C13  N^rKVEC 
C23  N«- ( \ N ) • . - \ N 

L 3 3 NC  iKVECl03«-“l 

L43  Z«-KVE£VN 

^DELCDDv 
v Z«-DEL  M 

C13  Z«-DELOP+.xM+.x*DELOP 

vMVNPR0BC03v 
v Z«-X  MVNPROB  Y » S f IND  f K 
C 1 3 -♦(  l = rrY)/CONTIN 

L"23  Y«~(0*,Y)/»Y 

C33  contin:s«-indel.  Y 
L43  IND«-  + /4Y 

C53  K>( ( (o2)*'l+rKVEC)xINDxx/Y)*~0.C 

C63  Z<-K'x*-0.5xX+.  xS+  ♦ xX 

v 

VINDELED3V 

v zmndel  x»sa;d;sc;sd 

C 1 3 ->(  l = rrX)/CONT 

C23  X*<0*,X)/>X 

C33  CONT  i SA«-KVEC/-rX 
C43  SAfSAo.xSrt 
L“53  D«-  + + /-fX 

C 6 3 SC<-f>AxD 

L'73  SD<-(\-l  + rX)«.-t~l  + rX 

C 8 3 SD<-SDxSA*0.5 

C 9 3 Z^SD-SC 

y 

vCONPROBCCn? 

? Z«-X  CONPROB  TM 
C 1 3 I«-l 

1 23  Z «-  ( rTM)C13rO 

C33  NEXrjZCI3^X  MVNPROB  TMCIJ3 

C43  I + «-l 

C53  -»(  IirZ)/NEXT 

v 

VUTPROBCD3V 
v Z«-PM  WTPROB  PV 
Cl  3 Z«-  + /PhxpV 

V 


A- 3 


